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1 Introduction 

In 1917 S. Kakeya posed the Kakeya needle problem: what is the smallest 
area which is required to rotate a unit line segment (a "needle") by 180 
degrees in the plane? Rotating around the midpoint requires ir/A units of 
area, whereas a "three-point U-turn" requires vr/S. In 1927 the problem was 
answered by A. Besicovitch, who gave the surprising answer that one could 
rotate a needle using arbitrarily small area. 

At first glance, Kakeya's problem and Besicovitch's resolution appear 
to be little more than mathematical curiosities. However, in the last three 
decades it has gradually been realized that this type of problem is connected 
to many other, seemingly unrelated problems in number theory, geometric 
combinatorics, arithmetic combinatorics, oscillatory integrals, and even the 
analysis of dispersive and wave equations. 

The purpose of this article is to discuss the interconnections between these 
fields, with an emphasis on the connection with oscillatory integrals and PDE. 
Two previous surveys ([^ and [Q) have focused on the connections between 
Kakeya-type problems with other problems in discrete combinatorics and 
number theory. 

These areas are very active, but despite much recent progress, our un- 
derstanding of the problems and their relationships to each other is far from 
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complete. Ideas from other fields may well be needed to make substantial 
new breakthroughs. 

2 Kakeya type problems 

Besicovitch's solution to the Kakeya needle problem relied on two observa- 
tions. The first observation, which is elementary, is that one can translate a 
needle to any location using arbitrarily small area. The second observation 
is that one can construct open subsets of of arbitrarily small area which 
contain a unit line segment in every direction. A typical way to construct 
such sets (not Besicovitch's original construction) is sketched in Figure 2; for 
a more detailed construction, see J^. 




Figure 1: To translate a needle, slide it by R units, rotate by roughly 1/-R, 
slide it back, and rotate back. This costs 0{1/R) units of area, where R is 
arbitrary. 

For any n > 2, define an Besicovitch set to be a subset of R" which con- 
tains a unit line segment in every direction. The construction of Besicovitch 
shows that such sets can have arbitrarily small measure in any dimension, 
and can even be made to be measure zero. Intuitively, this states that it is 
possible to compress a large number of non-parallel unit line segments into 
an arbitrarily small set. 

In applications one wishes to obtain more quantitative understanding of 
this compression effect by introducing a spatial discretization. For instance, 
one could replace unit line segments by 1 x 5 tubes for some < 5 <C 1 and 
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Figure 2: The iterative construction of a Besicovitch set. Each stage consists 
of the union of triangles. To pass to the next stage, the triangles are bisected 
and shifted together to decrease their area. 
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ask for the optimal compression of these tubes. Equivalently, one can ask for 
bounds of the volume of the 5- neighbourhood of a Besicovitch set. 

Rather surprisingly, these bounds are logarithmic in two dimensions. It 
is known that the (5-neighbourhood of a Besicovitch set in must have area 
at leastQ C/log(l/(5); this basically follows from the geometric observation 
that the area of the intersection of two 1x5 rectangles varies inversely with 
the angle between the long axes of the rectangles. Recently, U. Keich has 
shown that this bound is sharp. 

This observation can be rephrased in terms of the Minkowski dimension 
of the Besicovitch set. Recall that a bounded set E has Minkowski dimension 
a or less if and only if for every < 5 -C 1 and < £ -C 1, one can cover E 
by at most Ce5~°+^ balls of radius 6. From the previous discussion we thus 
see that Besicovitch sets in the plane must have Minkowski dimension 2. 

the analogous question in is known as the Kakeya conjecture: 
Does every Besicovitch set in R" have Minkowski dimension n? 

Equivalently, the Kakeya conjecture asserts that the volume of the 5- 
neighbourhood of a Besicovitch set in R" is bounded below by Cn,eS^ for any 
€ > and < 5 < 1. 

There are several variants of this conjecture (e.g. one could discuss Haus- 
dorff dimension instead of Minkowski dimension), but we shall not discuss 
these variants here for sake of exposition. 

The Kakeya conjecture is remarkably difficult. It remains open in three 
and higher dimensions, although rapid progress has been made in the last 
few years. The best known lower bounds for the Minkowski dimension at 
this time of writing is 

although I expect further improvements to follow very soon. 

One can discretize the conjecture. Let be a maximal (5-separated sub- 
set of the sphere S*""^ (so that Q has cardinality approximately and 
for each G f2 let T^^ be a 5 x 1 tube oriented in the direction u. The 
Kakeya conjecture then asserts logarithmic-type lower bounds on the quan- 
tity I [Ju^enTj- 

^Throughout this article, the letter C will denote various constants which vary from 
line to line. 
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The above formulation is reminiscent of existing results in combinatorics 
concerning the number of incidences between lines and points, although a 
formal connection cannot be made because the nature of the intersection of 
two S X 1 tubes depends on the angle between the tubes, whereas the in- 
tersection of two lines is a point regardless of what angle the lines make. 
However, it is plausible that one can use the ideas from combinatorial inci- 
dence geometry to obtain progress on this problem. For instance, it is fairly 
straightforward to show the Minkowski dimension of Besicovitch sets is at 
least (n+l)/2 purely by using the fact that given any two points that are a 
distance roughly 1 apart, there is essentially only one S x 1 tube which can 
contain them both. 

In the 1990s, work by J. Bourgain, T. Wolff, W. Schlag, A. Vargas, N. 
Katz, I. Laba, the author, and others pushed these ideas further. For in- 
stance, the lower bound of (n + 2)/2 for the Minkowski dimension was shown 
in 1995 by Wolff and relies on the 5-discretized version of the geometric state- 
ments that every non-dcgcncratc triangle lies in a unique two-dimensional 
plane, and every such plane contains only a one-parameter set of directions. 
However, there appears to be a limit to what can be achieved purely by ap- 
plying elementary incidence geometry facts and standard combinatorial tools 
(such as those from extremal graph theory). More sophisticated geometric 
analysis seems to reveal that a counterexample to the Kakeya conjecture, 
if it exists, must have certain rigid structural properties (for instance, the 
line segments through any given point should all lie in a hyperplane) . These 
type of ideas have led to a very small recent improvement in the Minkowski 
bound to (n + 2)/2 + 10~^°, but they are clearly insufficient to resolve the 
full conjecture. 

The Kakeya problem is a representative member of a much larger family 
of problems of a similar flavour (but with more technical formulations) . For 
instance, one can define a (3-set to be a subset of the plane which contains 
a /3-dimensional subset of a unit line segment in every direction. It is then 
an open problem to determine, for given f3, the smallest possible dimension 
of a /3-set. Low- dimensional examples of such sets arise in the work of H. 
Furstenburg, and it seems that one needs to understand these generalizations 
of Besicovitch sets in order to fully exploit the connection between Kakeya 
problems and oscillatory integrals, which we discuss below. Other variants 
include replacing line segments by circles or light rays, considering finite 
geometry analogues of these problems, or by replacing the quantity | U^^'^l 
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by the variant || Xli^ Xt^^ lip (the relevant conjecture here is known as the 
Kakeya maximal function conjecture) . Another interesting member of this 
family is the Falconer distance set conjecture, which asserts that whenever 
E is a. compact one- dimensional subset of R^, that the distance set {\x — y\ : 
x,y ^ E} is a one- dimensional subset of R. The discrete version of this is the 
Erdds distance problem - what is the least number of distances determined 
by n points? - and is also unsolved. For a thorough survey of most of these 
questions we refer the reader to 0; see also 0. 

3 The connection with arithmetic combina- 
torics. 

The Kakeya problem looks like very geometrical, and it is natural to ap- 
ply elementary incidence geometry facts to bear on this problem. Although 
this approach has had some success, it does not seem sufficient to solve the 
problem. 

In 1998 Bourgain introduced a new type of argument, based on arithmetic 
combinatorics (the combinatorics of sums and differences), which gave im- 
proved results on this problem, especially in high dimensions. The connection 
between Kakeya problems and the combinatorics of addition can already be 
seen by considering the analogy between line segments and arithmetic pro- 
gressions. (Indeed, the Kakeya conjecture can be reformulated in terms of 
arithmetic progressions, and this can be used to connect the Kakeya conjec- 
ture to several difficult conjectures in number theory such as the Montgomery 
conjectures for generic Dirichlet series. We will not discuss this connection 
here, but refer the reader to [|I|). 

Bourgain's argument relies on the following "three-slice" idea. Let VL and 

be as in the previous section. We may assume that the tubes are 
contained in a fixed ball. Suppose that | Uo^en-^^l comparable to 5" for 
some constant a; our objective is to give upper bounds on a, and eventually 
to show that a must be zero. 

By choosing an appropriate set of co-ordinates, one can ensure that the 
three shces 

Xt:={xeYC-^ ■.{x,t)e IJtJ 



6 



have measure comparable to 5" for t = 0, 1/2, 1. Because of the (5-discretized 
nature of the problem, one can also assume that the discrete set 




Figure 3: The left picture depicts six tubes pointing in different directions, 
and the three discretized slices Aq, A1/2, Ai, which in this case are three- 
element sets. The right picture depicts the set G associated to this collection 
of tubes. Note that the maps (a, b) — > a, (a, b) —>■ b, {a,b) —>■ a + b have small 
range (mapping to Aq, Ai, and 2^1/2 respectively), but the map (a, b) a—b 
is one-to-one. 

Morally speaking, every tube T,^ intersects each of the sets Aq, A1/2, Ai 
in exactly one point. Assuming this, we see that every tube T^^ is associated 
with an element of Aq x Ai. Because two points determine a line, these 
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elements are essentially disjoint as uj varies. Let G denote the set of all pairs 
of Aq X A\ obtained this way. Thus G has cardinality about 
The sum-set 

{a + 6 : (a, 6) G G} 

of G is essentially contained inside a dilate of the set Axj^'-, this reflects the 
fact that the intersection of with A1/2 is essentially the midpoint of the 
intersection of with A^ and A\. In particular, the sum-set is quite small, 
having cardinality only On the other hand, the difference set 

{a-h: {a,b) G G} 

of G is quite large, because the tubes T^^ all point in different directions. 
Indeed, this set has the same cardinality as G, i.e. about 

Thus, if a is non-zero, there is a large discrepancy in size between the 
sum set and difference set of G. In principle, this should lead to a bound on 
a, especially in view of standard inequalities relating sum sets and difference 
sets, such as 



A-B\ < 





B 
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\A\ 


\B\ 





(A summary of such inequalities can be found in [||). However, these argu- 
ments (which are mostly graph-theoretical) do not seem to adapt well to the 
Kakeya application, because we only are working with a subset G of x Ai 
rather than all of AqX Ai. 

To overcome this problem, Bourgain adapted a recent argument of T. 
Gowers which allows one to pass from arithmetic information on a subset of 
a Cartesian product to arithmetic information on a full Cartesian product. 
A typical result is: 

Theorem. Let A, B be finite subsets of a torsion free abelian group 
with cardinality at most N , and suppose that there exists a. set G G Ax B 
of cardinality at least aN"^ such that the sum set {a + h : (a, h) G G} has 
cardinality at most A^. Then there exists subsets A' , B' of A, B respectively 
such that A' — B' has cardinality at most a~^^N and A' , B' have cardinality 
at least a^N . 

Roughly speaking, this theorem states that if most of A + i? is contained 
in a small set, then by refining A and B slightly, one can make all oi A — 
B contained in a small set also. Such results are reminiscent of standard 
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combinatorial theorems concerning the size of sum and difference sets, but 
the innovation in Gowers' arguments is that the control on A' and B' is 
polynomial in a. (Previous combinatorial techniques gave bounds which 
were exponential or worse, which is not sufficient for Kakeya applications). 

Recently N. Katz and the author have obtained the bound (4n + 3)/7, 
by using control on both the sum-set Aq + Ai and the variant Aq + 2Ai, 
which corresponds to the slice A2/3. 

These results have remarkably elementary proofs. Apart from some ran- 
domization arguments, the proofs rely mainly on standard combinatorial 
tools such as the pigeonhole principle and Cauchy-Schwarz inequality, as 
well as basic arithmetic facts such as 

a + b = c + d a — d = c — b, 

a — b = {a — b') — {a — b') + (a' — b), and 

ao + 2bo = ai + 2bi, b'^ = b[^ai - b[ = 2(ao + 60) - 26i - (ao + 60). 

Further progress has been made by pursuing these methods, though it seems 
that we are still quite far from a full resolution of the Kakeya problem, and 
some new ideas are almost certainly needed. 

One possibility may be that one would have to use combinatorial esti- 
mates on product sets in addition to sum sets and difference sets, since one 
has control of {a + tb : (a, b) G G} for all t G [0, 1]. Discrete versions of such 
estimates exist; for instance, G. Elekes has recently shown the bound 

max(|A ■ A\, \A + A\) > Vl^^^ 

for all finite sets of integers A. However, these bounds do not adapt well 
to the continuous Kakeya setting because of the difficulty in discretizing 
both addition and multiplication simultaneously. A good test problem in 
this setting is the Erdds ring problem - determine whether there exists a 
(Borel) sub-ring of R with Hausdorff dimension exactly 1/2. This problem 
is known to be connected with the /3-set problem and the Falconer distance 
set problem. 

Interestingly, the Kakeya problem is also connected to another aspect 
of arithmetic combinatorics, namely that of locating arithmetic progressions 
in sparse sets. (A famous instance of this is an old conjecture of Erdos 
that the primes contain infinitely many arithmetic progressions of arbitrary 
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length). This difficulty arises in the Hausdorff dimension formulation of the 
Kakeya problem, and also in some more quantitative variants, because of the 
difficulty in selecting a "good" set of three slices in arithmetic progression 
in which to run the above argument. The combinatorial tools developed for 
that problem by W. Gowers and others may well have further applications 
to the Kakeya problem in the future. 

4 Applications to the Fourier transform 

Historically, the first applications of the Kakeya problem to analysis arose in 
the study of Fourier summation in the 1970s. 

If / is a test function on R", we can define the Fourier transform / by 

One then has the inversion formula 

Now suppose that / is a more general function, such as a function in the 
Lcbesgue space L^(R"). The Fourier inversion formula still holds true in the 
sense of distributions, but one is interested in more quantitative convergence 
statements. Specifically, we could ask whether the partial Fourier integrals 

SRfix) := / e2--V^(0 

converge to / in (say) norm (the pointwisc convergence question is also 
interesting, say for functions /, but this seems extremely difficult to show 
in two and higher dimensions. In one dimension this was proven in a fa- 
mous paper by L. Carleson). By the uniform boundedness principle, this is 
equivalent to asking whether the linear operators Sr are bounded on Z/'(R"') 
uniformly in R. By scale invariance it suffices to show this for Si: 

ll^i/ILnR") ^ ^II/IIl.(R")- 

The operator Si is known as the disk multiplier, because of the formula 

Sif = XbI 
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where B is the unit disk in R". In one dimension it is a classical result of 
Riesz that this operator is bounded on every L^, 1 < p < oo, and so Fourier 
integrals converge in LP norm. (Indeed, in one dimension Si is essentially 
the Hilbert transform). In higher dimensions, 5*1 is bounded in thanks to 
Plancherel's theorem, however the behaviour in is more subtle. One has 
an explicit kernel representation which roughly looks like 

/^±i\x-y\ 
(l + |^_^|)(n+l)/2 /(^) 

to be more precise, one must use Bessel functions instead of e='=*l^~^l. The 
kernel is only in LP when p > , so it was natural by duality arguments to 
conjecture that was bounded when < p < In 1971, however, C. 
Fefferman proved the surprising 

Theorem. If n > 1, then is unbounded on LP for every p ^2. 

In particular, one docs not have convergence for the Fourier inversion 
formula in higher dimensions unless p = 2. 

Roughly speaking, the idea is as follows. By duality it suffices to consider 
the case p > 2. Let i? be a large number, and let T be a cylindrical tube in 

with length R and radius y/R, and oriented in some direction ujt- Let 
ipT be a bump function adapted to the tube T, and let T be a shift of T by 
2R units in the cut direction. Then a computation shows that 

for all a; G T. 

To exploit this computation, one uses the Besicovitch construction to find 
a collection {T} of tubes as above, which are disjoint, but whose shifts T have 
significant overlap. More precisely, we assume that 

lU^i^Em 

T T 

for some K which grows in R (the standard construction gives K ~ log{R)/ loglog(i?). 
Then we consider the function 

/(a;) = 5]eTe^— -Vt(x), 

T 
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Figure 4: Four tubes Tj, their shifts Tj, and the wave packets ipi- The 
interference between the functions Siifji will cause the LP norm to be large 
when p> 2. 
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where er = ±1 are randomized signs. Using Khinchin's inequality (which 
roughly states that one has the formula | YIit ^tIt] ~ I/tP)^^^ with very 
high probability), one can eventually compute that 

Since K is unbounded, we thus see that 5*1 is unbounded. 

Feffcrman's theorem is an example of how a geometric construction can 
be used to show the unboundedness of various oscillatory integral operators. 
The point is while the action of these operators on general functions is rather 
complicated, their action on "wave packets" such as e^'^'^'^'^'^'4>T{x) is fairly 
easy to analyze. One can then generate a large class of functions to test the 
operator on by superimposing several wave packets together, and possibly 
randomizing the co-efficients to simplify the computation. 

The counter-example provided by the Besicovitch construction is very 
weak (only growing logarithmically in the scale Ft) , and can be eliminated if 
one mollifies the disk multiplier slightly. For instance, the counterexample 
does not prohibit the slightly smoother Bochner-Riesz operator S^, defined 

by 

smf = (1 - mxBioKo 

from being bounded for e > 0, because the analogous computation gives 

for all X e T. Indeed, the Bochner-Riesz conjecture asserts that Sf is indeed 
bounded on for all e > and 2n/{n + 1) < p < 2n/{n — 1). (For other 
values of p one needs £ > — || — |). This conjecture was proven by L. 
Carleson and P. Sjolin in 1972 in two dimensions, but the higher-dimensional 
problem is quite challenging, and only partial progress has been made so far. 
This conjecture would imply that the partial Fourier integrals will converge 
in if one uses a Ccsaro summation method (such as the Fejer summation 
method, which corresponds to e = 1). 

The Bochner-Riesz conjecture would be disproved if one could find a 
collection of disjoint tubes T for which the compression factor K had some 
power dependence on R as opposed to logarithmic, i.e. ii K > C~^R^ for 
some £ > 0. A more precise statement is known, namely that the failure of the 
Kakeya conjecture would imply the failure of the Bochner-Riesz conjecture. 
(More succinctly, Bochner-Riesz implies Kakeya). 
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In 1991, Bourgain introduced a method in which these types of imphca- 
tions could be reversed, so that progress on the Kakeya problem would (for 
instance) imply progress on the Bochner-Riesz conjecture. The key observa- 
tion is that every function can be decomposed into a linear combination of 
wave packets, by applying standard cutoffs both in physical space (by point- 
wise multiplication) and in frequency space (using the Fourier transform). 
After applying the Bochner-Riesz operator to the wave packets individually, 
one has to re-assemble the wave packets and obtain estimates for the sum. 
Kakeya estimates play an important role in this, since the wave packets are 
essentially supported on tubes; however this is not the full story since these 
packets also carry some oscillation, and one must develop tools to deal with 
the possible cancellation between wave packets. The known techniques to 
deal with this cancellation - mostly based on methods - are imperfect, so 
that even if one had a complete solution to the Kakeya conjecture, one could 
not then completely solve the Bochner-Riesz conjecture. Nevertheless, the 
best known results on Bochner-Riesz (e.g. in three dimensions the conjecture 
is known for p > 26/7 or p < 26/19) have been obtained by utilizing the 
best quantitative estimates of Kakeya type that are known to date. 

These techniques apply to a wide range of oscillatory integrals. A typi- 
cal question, the (adjoint) restriction problem, concerns Fourier transforms 
of measures. Let da be surface measure on (say) the unit sphere S"""^. 
The Fourier transform da of this measure can be computed explicitly using 
Bessel functions, and decays like Ixl"*-""^^/^ at infinity. In particular, it is in 
the class Lp(R") for all p > 2n/{n — 1). The restriction conjecture asserts 
that the same statement holds if da is replaced by fda for any bounded 
function / on the sphere. This question originally arose from studying the 
restriction phenomenon (that a Fourier transform of a rough function can be 
meaningfully restricted to a curved surface such as a sphere, but not to a 
flat surface like a hyperplane); it is also related to the question of obtaining 

estimates on eigenf unctions of the Laplacian on the torus (although the 
eigenfunction problem is far more difficult due to number theoretic issues), 
as well as obtaining estimates on solutions to dispersive PDE, as we shall 
see below. 

The restriction conjecture is logically implied by the Bochner-Riesz con- 
jecture, and is slightly easier to deal with technically. It has essentially 
the same amount of progress as Bochner-Riesz; for instance it is completely 
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solved in two dimensions, and is true for p > 26/7 in three dimensions. 
One uses the same techniques, namely wave packet decomposition of the 
initial function /, Kakeya information, and estimates to handle the can- 
cellation, in order to obtain these results. 

There are an endless set of permutations on these types of oscillatory in- 
tegral problems; more general phases and weights, square function and maxi- 
mal estimates, more exotic function spaces, bilinear and multilinear variants, 
etc. There are some additional rescaling arguments available in the bilinear 
case, as well as some L^-based estimates, but apart from this there are few 
effective tools known outside of Bourgain's wave packet analysis to attack 
these types of problems. 

One variant of the above problems comes from replacing Euclidean space 
by a curved manifold. There are some interesting three-dimensional examples 
of C. Sogge and W. Minicozzi which show that the Kakeya conjecture can fail 
on such manifolds, which then implies the corresponding failure of oscillatory 
integral conjectures such as the natural analogue of Bochner-Riesz. This may 
shed some light on the robustness of Kakeya estimates and their applications 
in variable co-efficient situations. Certainly the arithmetic and geometric 
techniques used to currently attack Kakeya problems do not adapt well to 
curved space. 

5 Applications to linear dispersive PDE 

We have seen in the previous section that Kakeya type problems have ap- 
plication to oscillatory integrals. One important type of oscillatory integral 
arises from the solutions to linear PDE. For instance, the solution to the free 
Schrodinger equation 

iut — Am = 0; m(0, x) = f\x) 
for X G R", t G R is given by the explicit formula 

for all t 7^ 0. We shall restrict our attention to the free (homogeneous) 
equation for sake of exposition, but much of this discussion can be extended 
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to the forced (inhomogeneous) problem 

iut — A.U — F{x, t); u{0, x) — 0, 



which is also of importance in applications. 

One way to see the connection between the free Schrodinger equation and 
the problems in the previous section is to take space-time Fourier transforms 



Prom these formulae one easily obtains that u must be the distribution 



with 6 being the Dirac delta. In particular, u is the inverse Fourier transform 
of a measure on the paraboloid r = 27r|(^p. 

Thus the problem of estimating the size of solutions to the Schrodinger 
equation is closely related to the restriction problem mentioned in the previ- 
ous section. By using the work of P. Tomas and E. Stein on the restriction 
problem, R. Strichartz showed a number of estimates for these types of equa- 
tions. A typical one is 



for all solutions to the two-dimensional Schrodinger equation with finite 
norm. Estimates of this type are known as Strichartz estimates; they en- 
capsulate certain smoothing properties of the solution (an function is 
smoother than an function) as well as decay properties in time (the spa- 
tial norm of u must eventually go to zero in order to make the space-time 
norm finite). They are essential in establishing the well-posedness and 
scattering theory of certain non-linear dispersive equations with low regular- 
ity initial data (see the next section), especially when combined with more 
elementary estimates such as conservation 




The Schrodinger equation then becomes 





«WIIl2(R") 



11/11 



for all t eR. 
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The above Strichartz estimate can be proven in a number of ways, either 
by a direct calculation using the Fourier transform, or by taking adjoints and 
using restriction theory, or by breaking u and / into wave packets and using 
some rudimentary Kakeya information. It has many generalizations; in fact, 
any linear equation which has a conserved energy and whose solutions decay 
in time will enjoy some family of Strichartz estimates. 

Strichartz estimates are one way of controlling the behaviour of linear 
PDE, but they are not the only such estimates. For instance, consider the 
problem of whether the explicit formula for the solution u{t,x) to the free 
Schrodinger equation actually converges back to the initial data f{x) when 
t — >■ 0. For test functions / this is easily verified, but the issue is more 
delicate when / is rough (cf. the discussion on the Fourier inversion formula 
above). A typical assumption is that / is only in a Sobolev space H'^ for 
some s G R, or in other words that 

UWhs := ||(1 + v^)7I|l^ 

is finite. (Roughly speaking, a function is in if the first s derivatives of 
the function are in L^. One can make sense of the operator (1 + \/— A)* by 
using the Fourier transform). 

One can ask the question of whether u{t, x) converges pointwise a.e. to 
f{x) as t ^ 0. This can be shown to be equivalent to estimating the maximal 
function 

sup \u{t, x) I 

i>0 

in some appropriate norm in terms of the norm of the initial data /. 
It is fairly easy to show that a.e. convergence fails for s < 1/4 and holds 
for s > 1/2, but the intermediate region requires a deeper understanding of 
the Schrodinger equation, and is not completely settled in two and higher 
dimensions. For instance, in two dimensions the best known positive result 
[^] is for s > 1/2 — 1/32. This argument relies on a bilinear estimate related 
to the restriction problem, which in turn is proven using wave packet de- 
compositions and Kakeya estimates. (One may view this type of argument 
as decomposing the Schrodinger wave into "particles" , and using Kakeya es- 
timates to control how often these particles collide with each other as time 
evolves) . 

The quantitative behaviour of the free Schrodinger equation is only par- 
tially understood at present. The Strichartz estimates are completely classi- 
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Figure 5: A schematic depiction of how a wave (such as the one drawn on 
the left) can be written as a superposition of wave packets or "photons". 
These objects are locahzed in space and also have a localized direction, and 
different wave packets are essentially orthogonal. There is no canonical way 
to perform this decomposition, but one usually uses a combination of spatial 
cutoffs and cutoffs in Fourier space. 
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fied, but these estimates do not give complete control on the solution. One 
basic issue (important in the study of energy concentration) is to under- 
stand the relationship between dispersion of u{x,t), and the distribution of 
the Fourier transform u(x,t) on the paraboloid {r = 27r|^p}. Heuristically, 
one expects u to disperse unless u{x, t) is concentrated in a small "cap" on 
the paraboloid, but obtaining a precise quantification of this phenomenon 
remains elusive, although estimates in this direction have been obtained and 
some interesting conjectures have been made. It seems likely that more pre- 
cise Kakeya estimates need to be developed before one has a satisfactory 
understanding of the size and dispersion of Schrodinger solutions. 

There is a close parallel between the theory of the Schrodinger equation 
and that of the free wave equation 

au{t,x)^0; u{0,x) ^ f{x),ut{0,x) ^0 

where □ = —d^ + A. To simphfy the exposition we restrict to the case of 
zero initial velocity. 

One can solve for u explicitly using the formula 

u{t) = cos{t^/^)f 

but this does not reveal much information about the size and distribution 
of u. The Fourier transform of u is supported on the hght cone {|r| = |^|}, 
which plays the same role as the paraboloid in the Schrodinger equation, 
although the cone has less curvature (so u has less dispersion), and the cone 
has bounded slope (so u has finite speed of propagation). 

A general class of problem is the following: given size and regularity 
conditions on the initial data / what type of size and regularity control does 
one obtain on the solution? 

From integration by parts (or from the above explicit formula) one has 
energy conservation 

J l\Vu{t,x)\' + l\u,{t,x)\' dx = j \\Vf{x)\' 

for all time t. This conservation law, and its generalizations, show that u has 
as much regularity as / when measured based spaces. However, control 
by itself does not reveal whether u focuses or disperses. To obtain better 
quantitative control on u one needs other estimates, such as W estimates. 
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The energy estimate is a fixed time estimate; it controls the solution at a 
specified time t. In the setting, fixed time estimates are very unfavorable, 
and require a lot of regularity for the initial data. A typical estimate is the 
decay estimate 

Ht)\\^^^^^^<c{i + \t\)-^--'y^ Y: I|vVL.(R") 

0<A;<s 

whenever s > (n + l)/2 is an integer. The necessity of this many derivatives 
is demonstrated by the focussing example, in which the initial data is spread 
out near a sphere of radius 1, and the solution u focuses (with an extremely 
high L°° norm) at the origin at time t = 1. 

However, one can obtain much better estimates, requiring far fewer deriva- 
tives, if one is willing to average locally in time. The intuitive explanation 
for this is that it is difficult for a wave to maintain a focus point (which 
would generate a large norm for p > 2) for any length of time. This 
phenomenon is known as local smoothing. Strichartz estimates are global 
smoothing estimates and thus fall into this category. A typical such estimate 
is 

in three spatial dimensions; this is the analogue of the Strichartz estimate 
for the two-dimensional Schrodinger equation mentioned earlier. Without 
the averaging in time, one would require 3/4 of a derivative on the right- 
hand side rather than 1/2; this can be seen from Sobolev embedding. These 
Strichartz estimates are usually proven by combining the energy and decay 
estimates with some orthogonality arguments. 

However, even Strichartz estimates lose some regularity. One may ask if 
there are estimates other than the energy estimate which do not lose any 
derivatives at all. Unfortunately, even if one localizes in time and assumes 
L°° control on the initial data, one still cannot do any better than control, 
as the following result of Wolff shows: 

Theorem. If n > 1 and p > 2, then the estimate 

II^IIlp([1,2]xR") ^ C'II/IIl°°(B{0,1)) 

cannot hold for all bounded / on the unit ball. 

The argument proceeds similarly to Fefferman's disk multiplier argument. 
Let {T} be a collection of disjoint tubes arranged using the Besicovitch set 
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construction as in Fefferman's argument, except that we rescale the tubes 
to have dimensions 1 x i?^^/^ rather than R x \fR. On each of these tubes 
T we place a "wave train" , which is basically e^^'^'^ times a bump function 
adapted to T. Let / be the sum of all these wave trains (although we may 
randomize the signs of these trains to simplify computations). 




I 

Figure 6: A schematic depiction of the Wolff example at time zero. Because of 
the zero initial velocity, the wave trains will move in two opposite directions. 

At time zero, the function / has low L°° norm. However, as time evolves, 
each wave train T splits as the superposition of two pulses, one moving in 
the direction ujt, and the other in the direction —cut- For times 1 < i < 2, a 
large portion of the wave train at T now lives in the shifted tube T. Because 
of the large overlap of these tubes, the norm of u is large for all 1 < t < 2; 
as with Fefferman's argument, it is about K^^^ p. By letting i? — >■ oo one 
can make K unbounded, and this gives the Theorem. 
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Figure 7: The Wolff example at a time 1 < t < 2; only the incoming wave- 
trains are shown. 

Because the Besicovitch construction has a logarithmic compression rate, 
one could get around this obstruction by requiring an epsilon of regularity 
on the initial data. The local smoothing conjecture of C. Sogge asserts that 
no further loss of regularity occurs, or more precisely that 

II«IIl.([1,2]xR") ^ CpM^ + V^)7Lp(R") 

for all e > 0,n(i ~ ^) " | and 2 < p < oo. This conjecture is easy when 
p = 2 or p = OO] the most interesting case is when p = 2n/ [n — 1). 

The local smoothing conjecture is extremely strong, and would imply 
many of the known estimates on the wave equation. It implies the Kakeya 
conjecture, for a counter-example to the Kakeya conjecture could be used 
to strengthen Wolff's argument to disprove the local smoothing conjecture. 
This conjecture also implies the Bochner-Riesz conjecture; the idea is to write 
the Bochner-Riesz multiplier Sf in terms of wave operators cos(tv— A) via 
a one-dimensional Fourier transform. However, the conjecture is far from 
settled; even in two dimensions, the conjecture is completely proven for only 
p > 74 (due to T. Wolff), and at the critical exponent p = 4 the conjecture 
is only known for e > 1/8 - 1/88 §, §. 
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There are several other wave equation estimates which are related to 
those discussed here. An active area of research is to obtain good bilinear es- 
timates on solutions to the wave equation, as opposed to the linear estimates 
described here; this will be discussed further in the next section. 

Since u can be written in terms of circular averages of /, there is a also 
close relationship between wave equation estimates and estimates for circular 
means. (Such circular means estimates can then be used, for instance, to 
make progress on the Falconer distance problem mentioned earlier). There is 
also an extremely strong square function estimate conjectured for the wave 
equation which, if true, would imply the local smoothing, Bochner-Riesz, 
restriction and Kakeya conjectures. It would also give estimates for other 
seemingly unrelated objects such as the helix convolution operator f ^ f * 
da, where da is arclength measure on the helix {(cost, sin t, t) : < t < 
27i} in R'^. (The connection arises because the Fourier transform of da is 
concentrated near the light cone). These estimates are quite difficult, and 
the partial progress which has been made on them has proceeded via Kakeya 
estimates. Although these deep wave equation estimates have not yet found 
significant applications, I am confident that they will do so in the near future. 

6 Applications to non-linear dispersive PDE 

Some of the estimates on linear PDE discussed in the previous section have 
been proven to be very useful in the analysis of non-linear perturbations of 
these PDE. We shall restrict our discussion mostly to non-linear wave equa- 
tions, although there has been parallel lines of research on similar equations 
such as non-linear Schrodinger and KdV. A detailed discussion of these topics 
can be found in the recent survey 

The idea that estimates on linear PDE lead to control of non-linear per- 
turbations is an old one, dating back at least to Picard. To illustrate the 
basic idea, consider the problem of solving the non-linear equation 

5x + e sin(x) = b, 

for X G R, where e, b are given, and e is small. Such an equation cannot be 
solved explicitly, however one can show that a solution exists, is unique, and 
depends analytically on the parameters e and b provided that e is sufficiently 
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small. To see this, simply rewrite the equation as 

X — 5~^6 + e5~^ sin(x). 

Since sin(a:) is Lipschitz, the right-hand side is a contraction if e is suf- 
ficiently small, and the claim then follows from the contraction mapping 
principle. More concretely, one constructs the solution x as the limit of the 
iteration 

^(n+i) 5-i6 + e5-^sin(a;(")) 

with x^^^ = (say). An essentially equivalent approach would be to expand 
x as a power series in e. 

This perturbative approach extends to non-linear evolution equations pro- 
vided that the non- linear equation is sufficiently well approximated by some 
linear or otherwise easily solvable equation (5a; = b, in the above example). 
This method is also very useful for demonstrating that qualitative phenom- 
ena associated to a special equation are robust under perturbations of that 
equation; for instance, one can use these techniques show that many non- 
linear wave equations obey the same type of properties (decay, resolution 
into plane waves, finite speed of propagation, energy conservation, few inter- 
actions between high and low frequencies, etc.) that the free wave equation 
enjoys. 

The iteration method gives existence, uniqueness, and analytic depen- 
dence on the data; these properties are collectively referred to as analytic 
well-posedness. In order to make the iteration converge, one often needs to 
localize in time, so that one only obtains local analytic well-posedncss. 

Iterative techniques do not extend well to strongly non-linear equations 
such as Navier- Stokes, but have been quite successful when applied to dis- 
persive equations such as the Korteweg-de Vries (KdV) equation and its 
variants, as well as non-linear Schrodinger and wave equations. The latter 
category includes such equations as the wave maps and Yang-Mills, as well 
as the more difficult quasi-linear equations such as the Einstein equations. 

A typical example is the semi-linear wave equation 

—Utt{t, x) + Au{t, x) = —u^{t, x);u{Q, x) = ef{x);ut{Q, x) = 

in three spatial dimensions, where we set the initial velocity to zero for sake 
of exposition, and the initial position to be small. This is one of the simplest 
non-linear perturbations of the free wave equation. 
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The ODE Uu = blows up in finite time, as one can see from such 
solutions as u{t) = \/2/{t — to). So one does not always expect a global 
solution, especially for large initial data. However, the Laplacian term in □ 
creates a dispersion effect, which can counteract this ODE blowup mecha- 
nism. Heuristically, if e is sufficiently small, then the dispersive effect of the 
Laplacian should overcome the blowup effect of the non-linearity, and one 
should obtain global existence. 

We can rewrite this equation in integral form as 

u{t) = £cos(iV^)/ + D-^u^{t) 
where 0~'^F is the unique solution u to the inhomogeneous wave equation 

Ou = F- m(0, x) = ut{0, x) = 0. 
An explicit formula for is given by Duhamel's formula: 

Consider the first few iterates of this equation. The first iterate would be 
the free solution uq := e cos{t\/—A)f. The second iterate would be obtained 
by inserting the free solution into the non-linearity, to obtain a better ap- 
proximation ui :— Uo + O~^ul. And so forth. It is clear that in order to make 
these approximations converge we must obtain estimates on such expressions 
as □~^-Uq. The norms are suited for this task, and by using the Strichartz 
estimates for the free wave equation (together with variants for the inhomo- 
geneous wave equation), one can indeed show analytic global well-posedness 
for data / in the Sobolev space H^^'^ for sufficiently small £, as the above 
integral equation turns out to be a contraction in a small ball in j. In fact, 
one can refine the analysis and also obtain scattering results (so that the non- 
linear solution eventually converges to a linear solution). The space H^^"^ is 
not chosen arbitrarily; it is the unique Sobolev space if* which is invariant 
under the natural scaling associated with the above equation, and is referred 
to as the critical regularity for this equation. Regularity smoother (resp. 
rougher) than critical is referred to as suhcritical (resp. supercritical) respec- 
tively. As a general rule, local well-posedness is relatively straightforward 
to establish for sub-critical regularities, quite delicate (and often unknown) 
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for critical regularities, and for super-critical regularities the techniques de- 
scribed above break down completely. For super-critical regularities one does 
not expect well-posedness, let alone analytic well-posedness, although it is 
quite plausible that solutions will exist and be unique for "generic" choices 
of supercritical data. 

There are several motivations for obtaining well-posedness results near 
the critical regularity, beyond the general benefit of testing, refining, and 
understanding the techniques used to prove these results. Firstly, the time of 
existence given by such results tends to grow as one approaches the critical 
regularity, and in the critical case one often gets global existence automati- 
cally for small data. Secondly, if blowup does indeed occur, the near-critical 
well-posedness theory often gives very precise control on the nature of the 
blowup (e.g. the energy must concentrate at a point). In certain cases, such 
control can then be combined with other facts (such as conservation laws) 
to ensure that blowup cannot actually occur. Finally, one can often refine 
the well-posedness analysis to obtain good control on various aspects of the 
solution, such as the transfer of energy from low frequency modes to high 
frequency modes or vice versa; this can also be used to prevent blowup and 
demonstrate global well-posedness even in cases where the energy is infinite. 

For semi-linear equations such as the one described above, a satisfactory 
theory has been constructed using the Strichartz estimates as the main tech- 
nical tool. The situation is slightly less satisfactory with equations which 
are still semi-linear, but have derivatives in the non-linearity. A typical 
example is the wave map equation for the sphere. Wave maps are maps 
: R""*"^ — > S"^ to a sphere which obey the non-linear wave equation 

□0 = -0(|V0r-|0,r). 

An important unsolved problem is whether wave maps with smooth data and 
finite energy remain smooth for all time in the two-dimensional case n — 2. 
(The result is false for n > 2 and true for n = 1, but the n = 2 case is far 
more delicate because the energy norm is critical in this case). 

If one tries to apply the above techniques, one is quickly faced with the 
need to estimate bilinear expressions such as 

Q(0,^) := V0 ■ - 0t ■ V^t. 

Most bilinear combinations of waves are strongest when 0, are both par- 
allel travelling waves, e.g. (p = cie**^^*^^^~*\ ip = 026**^^^^^"*^ However, the 
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above expression actually vanishes in this case. Because of this, the above 
expression is known null form. 

One could estimate this null form by the triangle and Holder inequalities, 
then using linear estimates such as the Strichartz estimates. However, this 
method does not take advantage of the cancellation in the null form and 
requires a fair amount of regularity on the solution (about half a derivative 
more than the critical regularity). In recent years it has been realized that 
in order to estimate these null forms accurately one needs genuinely bilinear 
estimates that take advantage of the null cancellation. 

As was observed already, null forms vanish when (p and ip have the parallel 
frequencies. Thus, one expects the largest contribution to a null form to come 
from interactions between waves which are transverse (i.e. their frequency 
vectors make an angle comparable to 1). The best tool we currently have to 
control such interactions are bilinear estimates. A typical such estimate 
is 

where 0, ip are free waves with initial position /, g respectively and ini- 
tial velocity zero, and Q is a null form. Strichartz methods can give this 
estimate, but with the right-hand side replaced by the symmetrized variant 
11/11^3/2 US' 1 1^3/2 5 which is less useful for applications. This particular estimate, 
together with some variants, was used by S. Klainerman and M. Machedon to 
demonstrate global well-poscdness of Yang-Mills and Maxwell-Klein-Gordon 
fields for finite energy data in three spatial dimensions. 

These estimates can be proven by taking Fourier transforms and using 
Plancherel's theorem. However, one can also prove these estimates using 
Kakeya methods. Roughly speaking, the idea is as follows. We shall induct 
on a scale parameter R; more precisely, we shall show that for each R > 
we have the bound 

for all space-time cubes of side-length R. Letting — > oo we obtain the 
result. 

Now fix R. We decompose the free wave (p into wave packets, each of 
which has a somewhat locahzed position and velocity. (This can be achieved 
by using cutoff functions in physical space to localize the initial data 0(0, x), 
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then using Fourier cutoffs to localize the direction of travel and the frequency 
magnitude). This decomposes (p into a superposition of functions '^rp(pT 
which are essentially supported on tubes near hght rays (see Figure 9) . (The 
size of these tubes depends on R and of the frequency magnitude of 0. For 
instance, if only contains frequencies of magnitude comparable to 1, then 
these tubes will have width about VR). 



t 




Figure 8: The wave packet decomposition. In this case consists of a single 
wave packet (pT- while ip is the superposition of two wave packets iJjt[ and 
ipT^. The parallel interaction between 4>t and 0^/ would be small because of 
the cancellation in the null form. 

We can decompose ip = J2t' similarly. We thus can split the origi- 
nal null form Q[(j),ijj) into a large number of smaller null form interactions 

(3(0T,^T')- 

One can divide these interactions into parallel interactions, in which T 
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and T' are parallel or nearly parallel, and transverse interactions^ in which 
T and T' intersect with angle comparable to 1. (If T and T' do not intersect 
at all, then Q{(t)T,'^T') is negligible). The exact dividing hne between what 
is parallel and what is transverse will depend on R and on the magnitude of 
the frequencies of (j) and ip. The parallel interactions are small because the 
null form vanishes when applied to parallel plane waves, so we concentrate on 
the transverse case. Now we use the basic geometric fact that two transverse 
tubes of length R and thickness ^fR can only intersect in a cube of size 
roughly \fR. To exploit this, we partition the original cube into smaller 
cubes q of side-length \/R, and observe that any pair (T, T') of tubes can 
only interact within a finite number of such smaller cubes. If one applies 
the induction hypothesis to each of q individually, and then sums up over g, 
one can (essentially) close the induction and recover the estimate at scale 
it! from the estimates at scale \fR. (One has to use some orthogonality 
properties of the 0t and iI^t' to properly do this). 

This induction on scales idea is mainly due to T. Wolff, who also extended 
the above argument to cover a near-optimal range of estimates in addition 
to estimates. These estimates and techniques are very recent, and have 
not yet found direct application to non-linear equations, but I am confident 
that they will do so in the near future. 

There are several active areas of research in non-linear wave equations. 
One important problem is to understand what happens at critical regularities. 
The iteration techniques described above seem to just barely fail to control 
the solution in many of these cases, and some new ideas are needed. Another 
problem is to extend the local existence theory given by iteration techniques 
to global existence. If there is a finite conserved quantity which controls the 
local time of existence, then this is straightforward, but a current area of 
research (by J. Bourgain, M. Keel, J. CoUiander, G. Staffilani, H. Takaoka, 
the author, and others) is to see what can be done if the conserved quantity 
is infinite. Another area, in which much progress has recently been made, 
is to extend the semi-linear theory to quasi-linear equations, in which the 
□ operator is replaced by a variable co-efficient operator depending on the 
solution (f). A typical equation is 

where gij is a smooth perturbation of the Kronecker delta 6ij, and is 
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Figure 9: The transverse interactions occuring at a single cube q of side- 
length y/R. By applying the induction hypothesis at this scale, and using 
orthogonality of wave packets to sum over all cubes q, one can close the 
induction. 
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some pseudo-differential operator of order -1. Tliis is a quasi-linear equa- 
tion of type similar to the Yang-Mills or Maxwell- Klein- Gordon equations. 
So far, progress on this problem (by H. Bahouri, J. Y. Chemin, D. Tataru 
and S. Klainerman) has come from adapting Strichartz estimates to variable 
co-efiicient situations when the metric is quite rough. A current area of re- 
search is to also adapt the bilinear estimates to the rough co-efiicient setting. 
For this the Fourier transform is no longer as efficient a tool, but it seems 
that physical space techniques such as Kakeya methods are more robust and 
should be able to handle this situation. 
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